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ABSTRACT 


Using a long (~150 ksec), broad-band (0.1-40 keV) BeppoSAX observation of 
the merging starburst system Arp 299 (=IC 694 -|- NGC 3690) we found the hrst 
unambiguous evidence of the presence of a deeply buried {Nh — 2.5 x 10^^ cm“^) 
AGN having an intrinsic luminosity of Tq.s-ioo keV — 1.9 x 10^^ erg s“^. The X- 
ray spectral properties of this AGN are discussed in detail as well as the thermal 
component detected at soft X-ray energies which, most likely, is associated with 
the starburst. 

Subject headings: galaxies: active - galaxies: individual (IG 694, NGG 3690) 

- galaxies: groups (Arp 299) - galaxies: nuclei - galaxies: starburst - X-rays: 
galaxies 


1. Introduction 

Studies of active objects at IR and X-ray wavelengths indicate that star-formation 
and AGN activity may be related (Fadda et al. 2002). The triggering mechanism for both 
phenomena could be the interaction or the merging of gas-rich galaxies. This generates 
fast compression of the available gas in the inner galactic regions, causing both the onset 
of a major starburst and the fueling of a central black hole raising the AGN activity (see 
Gombes, 2001 for a recent review). These two processes may proceed with different time 
scales and once initiated they probably have different lifetimes. Therefore a detailed study 
of the relative importance of starburst and AGN activity in the nearest (and brightest) 
objects may have profound impact on the understanding of the evolution and the fate of the 
gas in the interaction process, on the modeling the cosmological evolution of both AGNs and 
starburst objects (e.g. Franceschini, Braito and Fadda 2002), and on the evaluation of the 
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contribution of these sources to the energy density of the Universe (e.g. Fabian et ah 1998). 

The concomitant AGN and starburst activity is expected to happen in a high-density 
medium [Nh > cm“^), characterized by high dust extinction of the UV-optical 

flux and strong photoelectric absorption of the soft X-rays (e.g. Fabian et ah 1998). 
Therefore, the study of these active phases in galaxies becomes very difficult. Specihc 
examples are NGC 6240 and NGG 4945. Both objects have been classihed as LINER 
and/or starburst galaxies on the basis of optical (Veilleux et al. 1995) and mid-/far-IR 
(Genzel et al. 1998) spectroscopy. For both objects, however, the BeppoSAX PDS 
observations at E > 10 keV have clearly revealed the presence of a deeply buried AGN 
[Nh — few xlO^^ cm“^) with a QSO-like intrinsic luminosity in the case of NGG 6240 
(Lo. 5-100 keV ~ 6 X 10"^^ erg s“^; Vignati et al. 1999) and a Seyfert-like luminosity in the case 
of NGG 4945 (Lo. 5-100 keV ~ 7 x 10^^ erg s“^; Guainazzi et al. 2000). These two examples 
clearly show that optical and mid-/far-IR spectroscopy may not be sufficient to disentangle 
starburst activity from AGN activity, which is actually best probed in the hard [E > 6 keV, 
in order to sample also the Fe Ka line) X-ray energy band. 

To shed light on the starburst-AGN connection and its occurrence we have started a 
systematic and objective investigation in hard [E > 6 keV) X-rays of a flux-limited sample 
of IRAS galaxies. The sample consists of 28 galaxies selected from the IRAS Cataloged 
Galaxies and Quasars (see http;//irsa.ipac.caltech.edu/) as having /eo^m > 50 Jy or 
f 25 lira > 10 Jy. We stress here that no other selection criteria (e.g. established presence of 
an AGN, luminosities, IR colours, etc.) have been applied to the sample dehnition. 

In this letter we present and discuss BeppoSAX observations of Arp 299, a merging 
system (composed of IG 694 and NGG 3690, but simply quoted as NGG 3690 in the IRAS 
catalogue) located at D = 44Mpc (z = 0.011; Heckman et al. 1999 - H99), spectroscopically 
classihed as starbursting from optical (Goziol et al. 1998) and mid-/far-IR (Laurent et 
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al. 2000) observations. Its total FIR (43 — 123/im) luminosity, 2.86 x 10^^ Lq (following 
the recipe in Helou et al. 1988 and using the IRAS Faint Source Catalogue fluxes, 
Moshir et al. 1990), dominates the bolometric luminosity. As shown in this letter the 
BeppoSAX observations unveil a strongly absorbed {Nh — 2.5 x 10^^ cm“^) AGN with 
an intrinsic (unabsorbed) luminosity of Lo. 5-100 kev — 1-9 x 10^^ erg s“^: this is the 
first unambiguous evidence of the presence of an AGN in the interacting 
system Arp 299. For a detailed description of the system and for a summary of previous 
multiwavelength observations see Zezas, Georgantopoulos and Ward (1998; ZGW98), 
H99, Hibbard and Yun (1999), Alonso-Herrero et al. (2000), and Gharmandaris, Stacey 
and Gull (2002). In this paper we assume Hq = 75 km s“^ Mpc“\ and go = 0.5. 


2. Observations and data reduction 

Arp 299 was observed by BeppoSAX (Boella et al. 1997) in 2001 December 14 - 18 
for about 150 ksec. In this letter we use data collected from the Low Energy Goncentrator 
Spectrometer (LEGS), the Medium Energy Goncentrator Spectrometer (MEGS) and the 
Phoswich Detector System (PDS). 

The cleaned and linearized data produced by the BeppoSAX Scientihc Data Genter 
(SDG, see http://www.asdc.asi.it/bepposax) have been analyzed using standard software 
packages (XSELEGT vl.4, FTOOLS v4.2 and XSPEG vll.O). At the best spatial resolution 
of the BeppoSAX instruments (~ 2'), Arp 299 is not resolved; no signihcant source flux 
variation was detected over the whole observing period. In the PDS field of view (~ 60' 
radius) there are no known and bright (2-10 keV) X-ray sources except Arp 299. In the 
ROSAT All Sky Survey catalogue we found only two QSOs within 60' from Arp 299; given 
their lower flux and off-axis angle we can exclude them as sources of contamination in the 
PDS energy range. 
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To maximize statistics and S/N, the LEGS and MEGS source counts were extracted 
from a circular region of 4' radius; background counts were extracted from high-Galactic 
latitude “blank” fields (provided by the BeppoSAX SDG). The PDS spectrum extracted 
with the standard pipeline (with the rise-time correction applied) was provided directly by 
the BeppoSAX SDG; the simultaneously measured off-source background was used. In our 
analysis we took into account only data in the 0.1-4 keV and 1.8-10 keV energy range for the 
LEGS and MEGS respectively (as suggested in the BeppoSAX Gookbook, Fiore et ah 1999). 
Spectral channels corresponding to energies 10-40 keV have been used for the PDS data. 
LEGS and MEGS (PDS) source counts have been rebinned to have a S/N > 3 {S/N > 2) 
in each energy bin. Standard calibration hies released by the BeppoSAX SDG (September 
1997 version) have been used. The LEGS to MEGS and PDS to MEGS normalization 
factors were allowed to vary in the ranges proposed by the BeppoSAX Gookbook. Net 
exposure times (count rates) are: 63830 s (0.0087 ± 0.0004 cts/s) for the LEGS; 177050 s 
(0.0124 ± 0.0003 cts/s) for the MEGS; 75490 s (0.0852 ± 0.0187 cts/s) for the PDS. 

3. Spectral analysis 

Single-component models ^ do not provide an adequate description of the broad-band 
(0.1-40 keV) spectrum of Arp 299. A single-temperature thermal model or a single 
unabsorbed power-law model are both rejected at a conhdence level greater than 99.9%. 
Figure 1 shows the ratio between BeppoSAX LEGS, MEGS and PDS data and the best-ht 

^All models discussed here have been hltered through a foreground Galactic absorption 
of Nn^Gai = 9.92 x 10^® cm“^ (Dickey and Lockman 1990). The thermal component(s) have 
been modelled with the MEKAL model. All quoted errors are at the 90% conhdence level 
for 1 parameter of interest (Ay^ = 2.71). 
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unabsorbed power-law model (F ~ 1.9). A bump at ~ 0.8 keV, a line-like feature at 
~ 6.4 keV, and a big bump in the PDS energy range (10-40 keV) are clearly evident. The 
residuals suggest the presence of a soft thermal component around 0.8 keV, which is a 
characteristic signature in all known starburst galaxy X-ray spectra (e.g. Dahlem, Weaver 
and Heckman 1998). The simultaneous occurrence of a strong line-like structure 
at ~ 6.4 keV and of a big bump in the PDS energy range is the distinctive 
spectral signature of an obscured AGN (Matt et ah 2000). 

Given these features, we have performed a broad-band £t which includes: a) a soft 
thermal model; b) an unabsorbed power law with photon index F; cj a Gaussian emission 
line at ~ 6.4 keV, and d) an absorbed^ power law having the same photon index F. The 
combination of the spectral components b) and d) is usually called “leaky-absorber” model 
where the absorbed power law (accounting here for the PDS energy range) represents 
the “hrst order” AGN continuum emerging after transmission through an obscuring cold 
medium (torus? starburst related dust?). The unabsorbed power law (accounting here for 
the MEGS energy range) represents the primary AGN spectrum scattered into the line of 
sight by a warm, highly ionized gas located outside the absorbing medium. Furthermore, 
the residuals (not reported here) still show a clear line-like excess at ~ 3.5 keV, which 
has been modelled with a Gaussian emission line. The results of this best-fit model ^ 
(x^/d.o.f=104.2/118) are reported in Table 1 and are shown in Figure 2. 

^Given that the intrinsic is expected to be high (~ 10^^ cm“^) we have used the 
PLGABS model (Yaqoob 1997), which describes the X-ray transmission correctly taking 
into acconnt Gompton scattering. 

^A simpler model described by a power law plus a thermal component and a Gaus¬ 
sian line at ~ 6.4 keV is rejected by the data at a confidence level greater than 99% 
(xVd.o.f=162.1/122). 
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4. Discussion 
4.1. The AGN component 

The BeppoSAX data provide unambiguous evidence of the presence of AGN 
emission from Arp 299, with an intrinsic (i.e. unabsorbed) luminosity of Tq.s-ioo kev — 

1.9 X 10^^ erg s“^. 

The intrinsic power-law photon index (F = 1 . 79 ^ 003 ) is very similar to that 
usually observed in Seyfert 1 galaxies (Nandra et ah 1997). The primary continuum 
is highly absorbed {Nh = 2.5 x 10^^ cm“^, with a 99% conhdence level lower limit of 
~ 1.76 X 10^^ cm“^) and can escape from the absorber only at energies above 10 keV: 
this indicates that a deeply buried “Compton-thick AGN” lies in the interacting system 
Arp 299. The scattered flux fraction (given here as the ratio between the normalizations 
of the unabsorbed and absorbed power law components) implied from the hts is (5 ± 3)%, 
similar to that usually found in Seyfert 2 galaxies (e.g. Turner et ah 1997). 

At the spectral resolution of the BeppoSAX MEGS, the Gaussian line at ~ 6.4 keV is 
unresolved with a rest frame energy position of 6.42 ± 0.13 keV and an observed equivalent 
width (EW) of 636+PI; eV. This value is signihcantly larger than that observed in Seyfert 1 
galaxies (EW = 230±60eV, Nandra et ah 1997) but is similar to that measured in Seyfert 2 
galaxies (e.g. Bassani et al. 1999). The measured position is not consistent with He-like 
(6.7 keV) or H-like (6.96 keV) Fe, as expected if the line were produced inside the warm, 
highly ionized gas which scatters the primary AGN spectrum. However it is consistent with 
the low ionization Fe-Ka line from cold material. This in turn leads to two possibilities: 
the line could be produced by a reflected and/or a transmitted component^. The hrst 

^It should be recalled that Fe-Ka emission is also typical of High Mass X-ray Binaries 
(HMXBs) spectra (White et al. 1983). However, its equivalent width, either observed directly 



possibility however is unlikely since a significant refiection component cannot be easily 
accommodated within the BeppoSAX data (using the PEXRAV model, the 90% upper 
limit on the refiection fraction is 0.07). The second possibility invokes the production 
of the cold Fe-Ko; line by transmission through the same absorbing medium that affects 
the primary AGN continuum: indeed, the EW measured with respect to the transmitted 
component, ~ 7 keV, is consistent (within the errors on the intrinsic Nh) with what is 
expected from transmission (see e.g. Leahy and Creighton 1993), making this possibility 
highly plausible. We note here that the absorbing medium in composite starburst-AGN 
galaxies, such as Arp 299, is not univocally associated to the putative absorbing torus but 
the nuclear starburst itself could be a significant source of absorption (Fabian et al. 1998; 
Levenson, Weaver and Heckman 2001). The measured EW (with respect to the transmitted 
component) for Arp 299 line is a factor ~ 3-4 above the prediction from the model in 
Ghisellini et al. (1994) obtained using a simple toroidal geometry, which is a strong hint 
towards assuming the above argument apply. 

Finally, the unresolved Gaussian line at 3.381 q}3 keV (EW~ 125 eV) is consistent with 
the ArXVIII Ko; line. This line, which originates in the hot scattering medium itself (e.g. 
Netzer, Turner and George 1998), has been observed also in the Compton-thick Seyfert 2 
galaxy Circinus (EW~ 60 eV; Sambruna et al. 2001). Based on the flux ratio between the 
ArXVIII Ka line and the Fe He-like line (which should be produced inside the same hot 
scattering medium) observed in Circinus, in Arp 299 we would expect an Fe He-like line 
with EW ~ 400 eV, a value consistent (within a factor of 2) with the derived 90% upper 
limit. 


in HMXBs (typically EW ~ 0.3 keV: see White et al. 1983) or inferred for galaxies with 
X-ray emission dominated by HMXBs (Persic & Rephaeli 2002), is inconsistent with the 
large value observed here (if the emission has to be completely accounted for by HMXBs). 
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4.2. The Starburst Component 

The best-fit soft X-ray thermal component gives a kT ~ 0.86 keV. This value is 
consistent, within the errors, with that previously obtained by ZGW98 and H99 using 
ROSAT PSPC plus ASCA data. It is also similar to that usually found in other well studied 
starburst galaxies (e.g. Dahlem, Weaver and Heckman 1998). This thermal emission 
likely originates from the interaction between hot, low-density galactic winds and cold, 
high-density ISM (see H99). 

Starburst galaxies are also characterized by a hard 2-10 keV spectral component 
interpreted as being either thermal with kT ~ 5 — 10 keV or non-thermal with P ~ 1.5 — 2 
(e.g. Dahlem, Weaver and Heckman 1998). Such hard component is likely to be the 
integrated emission of X-ray binaries, and can be modelled by means of a cutoff power-law, 
f{E) oc (CPL model, see Persic and Rephaeli 2002). 

In order to check to what extent the 2-10 keV emission from Arp 299 can be related 
to the starburst, in our model we have replaced the scattered AGN component (labelled 
as b in section 3) with; a) a thermal component or b) the absorbed GPL model under 
the hypothesis that the main contribution is due to X-ray binaries. Both alternatives are 
statistically viable, with: a) the thermal model with kT = 6.5811:1^ or b) the unabsorbed 
GPL model with P = 1.65 ±0.15 {kT hxed at 8 keV). In both cases, the spectral parameters 
of the soft X-ray component and of the absorbed AGN component remain consistent 
(within the errors) with those reported in Table 1. Since the inferred 2-10 keV luminosity 
of both the thermal and the GPL components is ~ 2 x 10^^ erg s“^ (reasonable for a bright 
FIR galaxy like Arp 299 ) we cannot rule out that at least part of the observed 2-10 keV 
X-ray luminosity is due to the X-ray binaries directly associated with the starburst (which 
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may also contribute to the observed Fe-Ka emission) ^ . This conclusion was reached by 
ZGW98 and H99 using ASCA data. However the ASCA statistics and energy coverage were 
insufficient to clearly detect the 6.4 keV Fe-Ka emission line and to reveal the high energy 
bump The new broad-band BeppoSAX data instead clearly require the presence of an 
absorbed AGN component. 


5. Conclusions 

In a BeppoSAX observation of the merging starburst system Arp 299 the 0.1- 
40 keV data coverage clearly reveals, for the first time in this system, the presence 
of a deeply buried {Nh — 2.5 x 10^^ cm“^) AGN with an intrinsic (i.e. unabsorbed) 
luminosity of Lo. 5-100 keV — 1-9 x 10^^ erg s“F This AGN component was missed 
in previous multiwavelength observations. Assuming a Galactic standard value of 
Eb-v/AIh = 1-7 X 10“^^ mag cm^ (Bohlin et al. 1978) and the measured Nh value we argue 
that the AGN is completely absorbed in the optical and IR band {Ay > 1000 mag!!). 

The intrinsic AGN’s X-ray luminosity is a factor ~ 50 less than Lpin (~ 10^^ erg s“^). 
Thus the total FIR luminosity of the system cannot be entirely associated to the AGN 
even assuming an AGN UV luminosity a factor ~ 10 greater than the X-ray luminosity (as 
observed in QSOs, e.g. Elvis et al. 1994). This suggests that the bulk FIR emission may be 
due to the starburst, in agreement with the results obtained by Laurent et al. (2000) using 

^If this is the case, than the AGN scattered flux fraction reported in section 4.1, (5±3)%, 
should be considered as an upper limit. 

®We have directly verihed that the archival ASCA data are consistent with the best ht 
spectrum reported in Table 1 and are not deep enough to clearly detect the 6.4 keV Fe-Ka 
emission line. 
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5-16/im ISOCAM data. This conclusion may also be supported by the fact that the ratio of 
the soft X-ray thermal emission (likely associated with the starburst) to the FIR emission 
(~ 10“*^ in the case of Arp 299) is similar to other “bona hde” starburst galaxies such as 
M82 (~ 6 X 10-^) and NGC 253 (7 x lO"®). 

A question unsolved at the present stage is the location of the AGN inside the 
interacting system Arp 299. At the spatial resolution of the BeppoSAX instruments we 
are collecting photons from Arp 299 (= IC694 -|- NGC 3690) as well as from the nearby 
galaxies Arp296 and MCG-|-10-17-2a. It is safe to exclude the last two as the origin of 
the hard X-ray emission since they are very weak sources in both the radio and in the IR, 
and they were not detected in soft X-rays by ZGW98. On the other hand both IC694 and 
NGC 3690 could be the host of the AGN. We have recently retrieved the public Chandra 
and XMM-Newton data for this system. From a preliminary analysis a strong line at ~ 6.4 
keV is clearly present in NGC 3690 and maybe also in IC694, suggesting the possible 
presence of two AGNs. The results of this analysis, which will also include an estimate of 
the binaries contribution to the 2-10 keV emission and a detailed comparison of the X-ray 
emission with that at other wavelenghts, will be reported in a forthcoming paper. 

We thank V. Braito, P. Severgnini, T. Maccacaro, G. Trinchieri and the anonymous 
referees for useful comments. Partial hnancial support from ASI (I/R/073/01) and MIUR 
(Co£n-00-02-004) is acknowledged. 
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Figure Captions 


Fig. 1.— Ratio between the unabsorbed power-law model and the BeppoSAX LEGS (open 
triangles), MEGS (hlled circles) and PDS (open squares) data. The PDS to MEGS normal¬ 
ization factor is constrained to vary within the range indicated by the BeppoSAX Gookbook 
(Fiore et al. 1999). 

Fig. 2.— Data £t with a MEKAL thermal component, the “leaky-absorber” continuum and 
two narrow Gaussian lines (see Table 1). Panel a: LEGS (open triangles), MEGS (filled 
circles) and PDS (open squares) folded spectra and data; panel b: ratio between the data 
and the best-fit model; panel e: the unfolded model. Panel e insert: confidence contours 
(68%, 90% and 99% confidence level for two interesting parameters) for the power law photon 
index and the absorbing column density. 



Table 1: Results of the spectral fit (LECS+MECS+PDS): MEKAL thermal component + absorbed and scattered power 
law + narrow Gaussian lines (x^/d.o.f=104.2/118), and unabsorbed X-ray luminosities 




kT (keV)/r/EHne (keV) 

Norm 

Nh/EW (eV) 

Luminosity (erg 

s-i) 






0.5-2 keV 

2-10 keV 

10-100 keV 

SB 

MEKAL(“) 

0.86^;?^ 



1.12x10^^ 

7.41x1039 

5.36x1033 


Absorbed P. L. 

1 70 + 0-05 (c) 

' ^-0.03 

590l®23(rf) 





AGN 

Unabsorbed P. L. 

i.79iS:StM 



3.09x10^2 

4.93x10^2 

1.06x10^3 


Lines (rest frame) 

q 00 + 0.16 

n qq+0-30(g) 
U.Oy_o 26 

125/^® 






fi 42+0-13 
'0-++-0.13 

n n8+0-29(g) 

U.UO_o 26 

636^236 





NOTE: Errors are quoted at the 90% confidence level for 1 parameter of interest (Ay^ = 2.71). The LEGS to MEGS and PDS to MEGS 
normalization factors (0.7 and 0.95, respectively) are consistent with the known differences in the absolute calibration of the instruments. 
The total observed fluxes of Arp 299 are 8.16 x 10“^^ erg cm“^ s“^(0.5-2 keV), 1.13 x 10“^^ erg cm“^ s“^(2-10 keV) and 
3.23 X 10“^^ erg cm“^ s“^(10-100 keV). The metallicity of the thermal component was fixed to the Solar one. In units of 
[10“^®/(47ri2^) f neUH dV] @1 keV, where D is the distance of the source in cm, Ue and nn are the electron and H density in units of 
cm“^, and V is the volume filled by the X-ray emmitting gas in cm^. xhe two power law photon index have been fitted together. 

In units of 10“® photons keV“^ cm“^ s“^ @1 keV. Golumn density of neutral hydrogen (in units of 10^'^ cm“^) in addition 
to Nfi^Gai = 9-92 X 10^® cm“^. The lines are unresolved at the spectral resolution of the BeppoSAX MEGS. 

In units of 10“® photons keV“^ cm“^ s“^ at the energy of the line. 







Unabsorbed power low model 
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